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ABSTRACT 

Using a model for refractory clouds, a novel algorithm for handling them, and the latest gas- 
phase molecular opacities, we have produced a new series of L and T dwarf spectral and atmosphere 
models as a function of gravity and metallicity, spanning the Tog range from 2200 K to 700 K. The 
correspondence with observed spectra and infrared colors for early- and mid-L dwarfs and for mid- 
to late-T dwarfs is good. We find that the width in infrared color-magnitude diagrams of both the T 
and L dwarf branches is naturally explained by reasonable variations in gravity and, therefore, that 
gravity may be the "second parameter" of the L/T dwarf sequence. We investigate the dependence 
of theoretical dwarf spectra and color-magnitude diagrams upon various cloud properties, such as 
particle size and cloud spatial distribution. In the region of the L— >T transition, we find that no 
one cloud-particle-size and gravity combination can be made to fit all the observed data, though the 
role of binarity in the extant L and T dwarf data sets has yet to be fully determined. However, our 
results suggest that current ignorance of detailed cloud meteorology renders ambiguous the extraction 
of various physical quantities such as ToH' and gravity for mid-L to early-T dwarfs. Nevertheless, for 
decreasing Toff , we capture with some accuracy the transition from CO to CH4 in the K band, the 
FeH and CrH features at ~1.0 /imand ~0.85 /Ltm(at least for the early- to mid-Ls), the emergence of 
the neutral Na, K, Cs, and Rb alkali features, the growth of the Y/Z band peak, the disappearance 
of TiO and VO in L dwarfs, the evolution in the shapes of the J, H, and K bands, and the evolution 
of the K I features near 1.25 and 1.17 fj,m. We speculate that the subdwarf branch of the L 
dwarfs will be narrower in effective temperature, and, that for low enough metallicity, the L dwarfs 
will disappear altogether as a spectroscopic class. Furthermore, we note that the new, lower solar 
oxygen abundances of AUende-Prieto, Lambert, & Asplund (2002) produce better fits to brown dwarf 
data than do the older values. Finally, we discuss various issues in cloud physics and modeling and 
speculate on how a better correspondence between theory and observation in the problematic L— >T 
transition region might be achieved. 

Subject headings: general — stars: low-mass, brown dwarfs — radiative transfer — molecular processes — 
infrared: stars 



1. INTRODUCTION 

In the last ten years since the first identification of 
T and L dwarfs in the solar neighborhood, more than 
sixty T dwarfs and hundreds of L dwarfs have been 
classified spectroscopically and photometrically (Oppen- 
heimer et al. 1995; Delfosse et al. 1997; Burgasser et 
al. 1999,2000a,b,c; Kirkpatrick et al. 1999,2000; Geballe 
et al. 2002; Martin et al. 1999; Cruz et al. 2003; 
Knapp et al. 2004; Golimowski et al. 2004). Both in- 
frared and optical spectral features and colors have been 
used to define subtypes from LI to L9 and from TO to 
T9. These subtypes span an effective temperature (Toff ) 
range from ^2300 K to ^750 K and encompass the edge 
of the hydrogen-burning main sequence (Tefj ~ 1700 K; 
L5— L6). Water (H2O) bands dominate the spectra of 
both Ls and Ts. However, the L dwarfs are character- 
ized by the disappearance of the optical bands of TiO 
and VO (Jones & Tsuji 1997) that define the M dwarfs, 
the growth of near-infrared FeH and CrH features, the 
emergence of CH4 in the K (~2. 0—2.2 /zm ) and L' (~3.3 
fim) bands, and the progressive reddening of their J — K 
colors. The T dwarfs are distinguished by the appearance 
of strong CH4 features in the H (~1.5— 1.8 /zm) and the 
K bands, the weakening and disappearance of the metal 



hydride features, the growing importance of collision- 
induced absorption by H2 (Borysow & Frommhold 1990; 
Borysow, J0rgensen, & Zheng 1997), the strength of neu- 
tral alkali metal resonance lines shortward of 1.0 micron 
(Burrows, Marley, & Sharp 2000), and "anomalously" 
blue infrared colors. It is thought that the L/T transi- 
tion is near ~1300 K. All T dwarfs are brown dwarfs, 
with substellar masses, but many L dwarfs, however ex- 
otic, could be stars with masses above the main sequence 
edge. 

While most of these trends and characteristics are un- 
derstood in terms of general chemistry (Burrows & Sharp 
1999; Fegley & Lodders 1996; Lodders 1999), the central- 
ity of molecules in L and T dwarf atmospheres has put 
a new premium on the calculation of molecular opacities 
and abundances. Though theorists have gradually re- 
sponded to this challenge (Partridge & Schwenke 1997; 
Burrows et al. 2002a; Dulick et al. 2003; Burrows & 
Volobuyev 2003), their remain significant deficits in the 
relevant molecular databases. A glaring example is the 
absence of data on the hot bands of methane in the H 
band, but across the board molecular spectroscopy lags 
the atomic spectroscopy developed over the last hundred 
years to support classical stellar atmosphere calculations. 
As a result, theoretical models of L and T dwarf atmo- 
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the growing library of well-calibrated spectra, though 
possible in a crude sense, remains imprecise (Burrows et 
al. 2002b; Marley et al. 2002; Tsuji, Nakajima, & Yanag- 
isawa 2004; Burgasser, Burrows, & Kirkpatrick 2005). 

However, the most difficult challenge confronting the- 
orists designing models of L dwarfs and the L to T 
transition is the accurate incorporation of condensate 
clouds {e.g., calcium aluminates, silicates, iron). Such 
refractories form at the M to L transition temperatures 
(^2300— 1800 K), persist and thicken to dominate and 
define the L dwarfs, and when their atmospheric effects 
finally diminish are thereby responsible for the L to T 
transition (Burrows et al. 2001; Ackermann & Marley 
2001; Marley et al. 2002; Tsuji, Ohnaka, & Aoki 1999). 
Were it not for condensate clouds the Ls would not be 
red in infrared colors and get redder with decreasing Ted . 
In fact, the L dwarfs would not exist as a spectroscopic 
type ^. Hence, to understand the M— >L^T sequence and 
their spectra in detail requires a mastery of not only the 
chemical condensation sequences and the consequent ele- 
mental depletions ("rainout"; Burrows & Sharp 1999) in 
a gravitational field, but the spatial extent, particle size 
and shape distributions, grain optical properties, and me- 
teorology of clouds as well. These complications do not 
confront one who models most other types of "stars." 

After the first rush of basic discoveries, observers are 
now beginning to refine our knowledge of L and T dwarfs 
from the optical to the mid-infrared. Some target facts 
to explain include the brightening in the J (~1.2— 1.3 
/im) and Y/Z (~1. 0—1.1 fim) bands for the early T 
dwarfs (Dahn et al. 2002; Tinney, Burgasser, & Kirk- 
patrick 2003; Vrba et al. 2004), despite the decrease in 
bolometric luminosity along the spectroscopic sequence 
(Golimowski et al. 2004), the resurgence in the early Ts 
of FeH at ^^0.99 fim (McLean et al. 2003; Burgasser et 
al. 2004; Nakajima, Tsuji, & Yanagisawa 2004; Gushing, 
Rayner, & Vacca 2005), the non-equilibrium CO abun- 
dances inferred from spectral observations at 4.5—4.9 
/im (Golimowski et al. 2004), the appearance of the 2.2 
/tmand 3.3 /imCH4 features in the late Ls, and the nar- 
row apparent range of Toff at the M— >L and L— >T tran- 
sitions. In fact. Golimowski et al. (2004) suggest that 
Tcffis nearly constant from L7 to T4.5. Moreover, the 
actual prevalence of spectroscopic binaries among the ob- 
served L and T dwarfs, and the resulting misplacements 
on the HR diagram, remain a concern and a caution (Liu 
& Leggett 2005; Gefino & Kulkarni 2005; Burgasser et 
al. 2005). 

The brightening in the J band is the most intrigu- 
ing anomaly. Its explanation must be the rapid thin- 
ning out of the clouds in the spectrum-forming region 
of the brown dwarf atmosphere during the L— ^T transi- 
tion. But how the effective opacity of the silicate clouds 
decreases so quickly with T^s and spectroscopic subtype 
to yield heavy-element depleted T dwarf atmospheres has 
yet to be explained. Ideas include the emergence of holes 
in the clouds, whose areal filling factor would need to in- 
crease with spectral subtype (Burgasser et al. 2004), 
a sudden collapse of the cloud deck (Tsuji & Nakajima 
2003; Tsuji, Nakajima, & Yanagisawa 2004), or an in- 

^ This fact serves to emphasize that the span and importance of 
the L dwarf family is a function of metaUicity, and for subdwarfs 



crease in the "sedimentation efficiency" of the clouds, 
with a concommitent rapid increase in the silicate parti- 
cle size (Knapp et al. 2004) . All these models have their 
strengths and weaknesses, and all are ad hoc, with lit- 
tle physical or meteorological justification for the rapid 
transitions proffered. 

There is emerging evidence that the L and T spec- 
troscopic sequences are not a one-parameter family in 
Toff , but that there is a second (and, perhaps, a third) 
parameter. Burrows et al. (2002b) and Knapp et al. 
(2004) suggest that this second parameter is gravity and 
our calculations here lend strong support to this notion. 
The widths of both the L and the T dwarf regions in 
HR diagrams are best explained with a range of gravi- 
ties. MetaUicity may also play a role, but a subordinate 
one, except for rare outliers, and we explore in this pa- 
per the metaUicity dependence of L and T dwarf model 
atmospheres, spectra, and colors. 

In this paper, we provide a new generation of L dwarf 
and T dwarf model spectra, colors, and atmospheres. We 
use updated FeH, GrH, TiO, and VO gas-phase opacities 
and abundances and investigate the dependence on cloud 
model parameters, in particular particle size and cloud 
extent/shape. To do so, we have developed a sophisti- 
cated and self-consistent algorithm for handling clouds 
in atmosphere calculations. With this algorithm, we 
demonstrate the ambiguities that remain in any treat- 
ment of the L— >T transition, while generating good spec- 
tral and color models for both the early- and mid-L 
dwarfs and the mid- and late-T dwarfs. In the process, 
we provide and discuss new spectral models, temper- 
ature/pressure profiles, color- magnitude diagrams, and 
color-color diagrams for a broad range of Toff s , gravi- 
ties, metallicities, and particle sizes. 

In 321 we list the numerical tools and databases we have 
used to generate self-consistent molecular atmospheres 
and discuss the numerical method we employ to handle 
clouds and convection in L and T dwarf atmospheres. 
We review in ^ general cloud physics and our cloud 
model parametrizations. Then, we continue in 2] with 
a discussion of the effects of cloud shape, particle size, 
and gravity. The associated models are provided only 
to explore parameter dependences; they are not our best 
model results. We then introduce our baseline model, 
which, though manifestly limited given the current rudi- 
mentary knowledge of cloud physics and incorporating 
as it does simple ansatze (as do all other published mod- 
els), nevertheless gives acceptable fits to the observations 
outside of the problematic L9-T3 spectral range. This 
model suite is being provided to astronomers to aide 
them in understanding the rich harvest of L and T dwarf 
spectral data. In [J^lwe discuss representative tempera- 
ture/pressure (T/P) profiles that serve as the jumping- 
off point for further speculation on the L^T transition. 
Next, in iJHlwe provide spectra for our fiducial model set 
as a function of Toff , gravity, and metaUicity, as well as 
a few comparisons between our spectral models and data 
to establish that the models fit reasonably well. This is 
followed in iJZlby Mj versus J — K, Mk versus ,/ — K, 
and Mz' versus i' — z' color-magnitude diagrams for these 
baseline models. For comparison, data for M, L, and T 
dwarfs with known parallaxes are superposed. Finally, 
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2. NUMERICAL TOOLS AND CHEMICAL DATABASES 

In a study such as the one we conduct for this paper 
many self-consistent spectra and atmospheres must be 
generated for a variety of physical parameters. To do 
so expeditiously, it is undesirable and expensive to cal- 
culate "on the fly" the opacities at a given wavelength 
from the extensive line lists available. Rather, we pre- 
compute and tabulate at 5000 frequency points, logarith- 
mically spaced from 0.3 /xmto 300 yum, the absorption 
and scattering opacities on a mesh in T/P space for ^30 
of the dominant molecular and atomic species found in 
L and T dwarf atmospheres. These tables are then used 
with a table of the corresponding molecular and atomic 
abundances at the desired metallicity to produce a large 
table of total opacities as a function of frequency, tem- 
perature, and pressure. Such a table fully incorporates 
the physics of line broadening and stimulated emission 
and is interpolated in during the iteration of the spec- 
tral/atmosphere model until a converged T / P profile and 
a spectrum in radiative equilibrium are obtained. In- 
stead of requiring hours on a supercomputer to generate 
a model, this procedure typically requires only minutes 
on a single-processor workstation. To obtain spectra at 
higher frequency resolution, the derived T/P profile can 
be reused with a higher-resolution opacity table, without 
the need to iterate further. 

The atmospheric structure is computed using a spe- 
cific variant of the stellar atmospheres code TLUSTY 
(Hubeny 1988; Hubeny 1992; Hubeny & Lanz 1995), 
with modifications as described in Burrows et al. 2002b, 
Sudarsky, Burrows, & Hubeny (2003) and Hubeny, Bur- 
rows, & Sudarsky (2003). We have developed a compu- 
tational scheme to treat clouds and their interplay with 
convection zones efficiently. The approach will be de- 
scribed in detail elsewhere (Hubeny et al., in prep.); here 
we briefly mention some of its most important aspects. 

We have reformulated the general linearization scheme 
in the so-called Rybicki-type approach (Mihalas 1978), 
which consists in writing the global linearization ( Jacobi) 
matrix in block form, where each block contains an indi- 
vidual state parameter for all depth points. This allows 
us to avoid using Accelerated Lambda Iteration (ALI) to 
treat the radiative transfer equation. Instead, we use full 
linearization of the radiation intensities, which leads to 
much faster and more stable convergence, without com- 
promising on computer time. 

Convection is treated in the mixing-length formalism, 
with a mixing length of one pressure scale height. In 
actual brown dwarf atmospheres, the temperature gra- 
dient is very close to the adiabatic gradient. During the 
linearization iteration, it often happens that the current 
gradient is somewhat lower than the adiabatic gradient, 
even when in reality it should be larger. We have devised 
a procedure that detects such spurious "convection gaps" 
and recomputes the temperature structure in such a way 
that at those points the convective flux is indeed non-zero 
before entering the next iteration of the linearization. 

In addition to introducing a flexible way to treat cloud 
shapes (see ti3.1|l , we have developed several ways to self- 
consistently determine the cloud's position. It turns out 
that it is more stable not to recalculate the cloud posi- 
tion after each iteration; instead it is preferrable to per- 



the clouds fixed between updates. In addition, we have 
introduced a flexible depth-point rezoning scheme. Es- 
sentially, we place additional depth points close to the 
cloud base whenever the new cloud position is recalcu- 
lated. Finally, in some cases the cloud position is found 
to oscillate between two (sometimes even three) distinct 
positions. This problem is avoided by setting the up- 
dated position of the cloud base between the previous 
position of the cloud base and the intersection of the T/P 
profile and the condensation curve. At convergence, the 
cloud base will be at the proper intersection point. 

The equation of state from which we obtain the 
T/P/density(p) relation is that of Saumon, Chabrier, 
& Van Horn (1995). The molecular compositions and 
ionization fractions (Burrows, Sudarsky, & Lunine 2003) 
are obtained using the chemical equilibrium code SOL- 
GASMIX (Burrows and Sharp 1999), with updated ther- 
mochemical data for Ti and V compounds, silicates and 
calcium aluminates, H~, and metal hydrides. The abun- 
dance code incorporates prescriptions for rainout and de- 
pletion due to condensate formation of refractory sili- 
cates, aluminates, titanates, iron, water, and ammonia. 
Despite the overall complexity of the problem, the most 
important atoms and non-refractory molecules comprise 
a small subset: H2, H2O, CH4, CO, N2, and NH3, FeH, 
CrH, TiO, VO, Na, and K. 

Our always-evolving opacify/spectroscopy database 
would require a paper in itself to explain, but is par- 
tially described (with references) in Burrows et al. 
(1997,2001, 2002a,2003). The opacities of the alkah metal 
atoms are taken from Burrows, Marley, & Sharp (2000), 
which are similar in the line cores and near wings to 
those found in Burrows & Volobuyev (2003). For future 
work, we have extended the opacity database into the 
UV down to ~0.08 /^m. Our solar- metallicity elemen- 
tal abundances are taken from AUende-Prieto, Lambert, 
& Asplund (2002), which supercedes Anders & Grevesse 
(1989). AUende-Prieto, Lambert, & Asplund have found 
that in the past the solar abundances of oxygen and car- 
bon were overestimated by as much as ^35—45%. While 
this change has but subtle consequences for stellar and 
solar atmospheres, since the water abundance so dom- 
inant in L and T dwarf atmospheres scales almost di- 
rectly with the oxygen abundance, these atmospheres 
are much more directly affected by the change. Using 
Spitzer IRAC data and our theoretical models, Patten 
et al. (2004) show that the new abundances fit the mid- 
infrared data much better. Hence, and in a curious and 
unexpected way, the new solar abundances are partially 
validated by brown dwarf observations. 

Cloud physics, structure, character, and opacities play 
a major role in the M— i>L transition, L dwarf atmo- 
spheres, and the L^T transition and so are discussed 
separately in the following section. 

3. CLOUDS AND CLOUD MODEL DESCRIPTION 

In the atmospheres of late M dwarfs with Tcff ^2500 
K, the most refractory compounds start to condense. 
These are mostly the calcium aluminum oxides grossite 
(CaAUOr = CaO + 2(Al203)) and hibonitc (CaAliiOig 
= CaO + 6(Al203), among others, and first appear in 
the upper radiative, not the lower convective, zones. Fig- 
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dwarfs, L dwarfs, and brown dwarf atmospheres. As Fig. 
^indicates, this first phase of condensate cloud formation 
is foUowed, as the temperature drops, by the condensa- 
tion of titanium compounds, at the expense of gas-phase 
TiO, and later by the condensation of VO. Therefore, the 
disappearance of TiO and VO lags slightly, but follows 
close on the heals of, the first emergence of refractory 
clouds. Due to the low elemental abundances by num- 
ber of titanium ('^8x10"^) and vanadium (~10^*), the 
compounds into which TiO and VO condense (see Fig. 
^ are themselves not major components of this radiative 
cloud layer. 

In retrospect, this early transition from late M dwarfs 
to L dwarfs was first identified by Jones & Tsuji (1997), 
who for late Ms noted the near simultaneous weakening 
of the TiO and VO spectral features with the shallowing 
of the water troughs. The latter is naturally explained 
by the appearance of a new continuum absorber, clouds 
of the refractory calcium aluminates. However, though 
oxygen is abundant, the elemental abundances of calcium 
and aluminum are not, just ~2.0xl0-6 and ~2.8xl0"^ 
respectively. As a result, this first generation of clouds is 
rather thin, though the particle sizes in radiative zones 
can be small (^0.5— 5 microns. Cooper et al. 2003). It is 
not until the appearance at slightly lower temperatures of 
condensates containing silicon, magnesium, and/or iron, 
with elemental abundances '^10 times those of calcium 
and aluminum, that the areal mass density of refractory 
clouds can be respectable. With the lowering of Tcff and 
the appearance of condensed iron and silicates (such as 
forsterite, Mg2Si04, and enstatite, MgSiOa), the photo- 
spheric layer begins to coincide both with the thickening 
silicate cloud layer and with a convective zone, and the 
pattern of the mid-L dwarf atmosphere is firmly estab- 
lished. Therefore, the transitions from M to L dwarf, 
from calcium aluminate to silicate/iron condensates, and 
from radiative to convective clouds overlap, though not 
perfectly. 

As Fig. ^ indicates, at a given pressure the many 
condensates can appear in a narrow range of tempera- 
ture. Furthermore, whatever the particle sizes and opti- 
cal constants, the optical depths of such clouds are them- 
selves sufficient after the early Ls to trip convection (and 
the associated updrafts and downdrafts) where there are 
clouds. Hence, after the early Ls (for which the first con- 
densates inhabit a stably-stratified radiative zone) ev- 
ery condensate whose condensation curve intersects the 
dwarf's T/P profile will most probably reside in a com- 
mon convection zone. It does not make sense to assume 
that each condensate is a separate, isolated layer, like a 
"pousse cafe." Rather, as the Tos decreases and the first 
clouds thicken, tripping convection in the atmosphere, 
condensation and grain growth will occur over the same 
finite time and pressure/radius range for many condens- 
ables. The kinetics of such a soup of growing condensates 
poses a daunting problem never before addressed, and 
not addressed here, though grain kinetics in the brown 
dwarf context has been receiving some attention of late 
(Helhng et al. 2001, 2004; Woitke & HeUing 2003,2004). 
Nevertheless, the equilibrium particle size distribution, 
achieved through the balance of growth processes in the 
convective zone and grain destruction at and below the 



T/P-profile/condensation-curve intercept, is poorly con- 
strained by theory (Cooper et al. 2003; Ackerman & 
Marley 2001). Furthermore, the optical constants of het- 
erogeneous grains of indeterminate composition and lay- 
ering are not easily derived from first principles. Qualita- 
tively, it is clear that the modal particle size of grains in 
stable radiative zones is smaller (~0. 1—5.0 ^m) than in 
turbulent convective zones (^10—150 /im) (e.g., Cooper 
et al. 2003), but confidence in the current analytic esti- 
mates should not be great. 

Once formed, clouds will not extend throughout the 
atmosphere above (to lower pressures and temperatures) 
the condensation intercept, but will settle into a tighter 
spatial distribution bounded (approximately) from be- 
low by this intercept. Hence, the upper atmosphere will 
be depleted of the heavy elements (such as Si, Mg, Fe, 
Ca, and Al) of which condensates are comprised. The 
T dwarfs themselves clearly indicate the general cor- 
rectness of this conclusion; their clear spectra demand 
that their condensate clouds reside at great depths be- 
low the formation region of the emergent spectrum (Mar- 
ley et al. 1996; Burrows et al 2002b). The same is 
true of Jupiter and Saturn, whose atmospheres are de- 
void of heavy elements (Fegley & Fodders 1996), but 
which must house at depth (^1000—3000 bars) both sil- 
icate and iron clouds. The oxygen fraction in an L or T 
dwarf is so large (~5xl0~^, solar) that condensation de- 
pletes its upper-atmospheric abundance by no more than 
^^15% (Burrows & Sharp 1999). Below the cloud base 
(at higher pressures) is the "infinite" reservoir of heavy 
elements that extends throughout the dwarf and which 
sets the heavy-element abundance boundary condition at 
the cloud base. However, the heavy elements that may 
have once existed in the upper atmosphere before con- 
densation do not all remain in the cloud once formed. 
How much cloud material does remain in the cloud de- 
pends on the dynamics of the cloud itself. The actual 
areal mass density, spatial extent, and shape of a cloud 
is not easily determined and depends upon the interplay 
of turbulent gas motions, grain growth, overshoot, tur- 
bulent and eddy diffusion, and gravitational settling, as 
does the grain size distribution itself. Ackerman & Mar- 
ley (2001) and Cooper et al. (2003) suggest that the 
thickness of a cloud scales with the local pressure scale 
height and that the cloud material is distributed expo- 
nentially from its base with a scale height between ~0.3 
and ~1.0 pressure scale heights. This seems reasonable, 
but the effect of multiple condensation intercepts for the 
different species and of an extended convection zone that 
can encompass these multiple condensation lines is not 
at all clear. Nor is it clear how the modal particle size 
varies with altitude, though there are some hints (Cooper 
et al. 2003). 

For the time being, this state of complexity demands 
simple parametrizations. Our approach has been to 
parametrize the spatial extent and shape of the conden- 
sate cloud and to test various composite optical proper- 
ties and modal particle sizes, loosely guided by extant 
theory and previous practice (Cooper et al. 2003; Bur- 
rows et al. 1997,2001; Ackerman & Marley 2001; Al- 
lard et al. 2001) ^. For a given parametrization, we 

* For a eiven modal Darticle size, we use the Deirmendiian (1964) 
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have compared the resulting color-magnitude diagrams 
and spectra with observations in an attempt to settle on 
the fiducial condensate-cloud approach that can be used 
for our next generation of atmosphere/spectral models 
(pubhshed herein) for TeffS from 2200 K to 700 K and a 
variety of gravities and metallicities, collectively encom- 
passing the L and T dwarf ranges (This is our model E 
described in H3.1|l . Since the mid- to late-Ts are not af- 
fected by such clouds, this effort is germane only to L 
dwarf models and the L— s-T transition. In the process, 
we have determined the dependence of dwarf spectra and 
colors on particle size and cloud morphology. 

3.1. Cloud Parametrizations 

We assume that the clouds are uniformly distributed 
over the sphere and that there are no holes in either lati- 
tude or longitude. The position of the high-pressure base 
of the cloud is set at the intersection of the T/P profile 
and the condensation curve. The corresponding quan- 
tities are denoted by Pq and Tq, and the corresponding 
column mass by rric- (The pressure at a given depth is 
equal to g x to, where g is the surface gravity and the 
areal mass column density is to.) The opacity (per gram 
of atmospheric material) of a condensate a is given by 

Ka{y, m) = Na Ma — Sa Cra{v, flo) /(to) , (1) 

M 

where Na is the mixing ratio of species a. Ma is its molec- 
ular weight (in grams), ^ is the mean molecular weight 
of the atmospheric material, Ajv is Avogadro's number, 
Sa is the supersaturation ratio, aa(i^, ag) is the opacity 
per gram of species a for frequency v and modal particle 
size/radius ao, and /(to) is our cloud shape function, de- 
fined such that the full presence of a cloud is denoted by 
/ = 1 while its complete absence is denoted by / = 0. 

Na and Ma are taken as known quantities (depending 
on the assumed composition and varying with metallic- 
ity). The supersaturation parameter and the modal par- 
ticle size are taken as input parameters of the model; for 
iron Sa is 0.01 and for the silicates such as forsterite it 
is 1.0. The cross sections for both true absorption and 
scattering are taken from precalculated Mie tables, in- 
terpolating to the assumed particle size oq. 

In the absence of an ab-initio theory of cloud forma- 
tion and merger, we envisage a simple empirical model 
of a composite cloud having a flat part just below the 
intersection point (upper cloud deck), and having an ex- 
ponential decline on the both sides of this flat part. The 
extent of the flat part, and the exponents of the decline 
on both sides are free parameters of the problem. The 
cloud shape function is thus given by 

/(to) ~ [m/miY^ , to < nil , 

/(to) — 1 , TOi < m < TOq , 

/(to) = (to/toq)"'^ , TO > Too , (2) 

where toi and toq are the upper (low-pressure) and lower 
(high-pressure) boundaries of the flat part of the cloud. 
si and S2 are indices such that 1/si and l/s2 times the 
local pressure scale heights are the scale heights of the 
particle number distributions in the exponential tails. 
TOi can be parameterized in different ways. For a sin- 



pressure scale height would indicate the cloud's approx- 
imate extent. However, for multiple cloud condensates 
with a range of TqS that inhabit a convectively mixed 
region, the flat portion should encompass the multiple 
cloud T/P intercept points, with an extent roughly in- 
dicated by the thickness of the condensate band in Fig. 
n This is why we have introduced the flat portion in eq. 
0. 

To study the influence of the cloud shape function on 
the temperature structure and the predicted spectra, we 
have considered the following five cases. In model A, the 
cloud declines exponentially below the high-temperature 
intersection point, but extends upward all the way to the 
surface (si = 0; S2 = 10). Model B corresponds to having 
the cloud everywhere in the atmosphere (si = 0; S2 = 0). 
Models A and B are similar and roughly correspond to 
the "Dusty" models of Chabrier et al. (2000) and AUard 
et al. (2001). 

Model C corresponds to a cloud decaying above and 
below the intersection point without having a flat part 
(si = 6; S2 = 10); model D represent a cloud that decays 
above the intersection point (si = 6), while extending all 
the way down deep into the dwarf (s2 = 0); and, finally, 
model E corresponds to our generic, fiducial cloud - a 
cloud with a flat part between the intersection point of 
one of the most refractory condensates (such as grossite) 
and the intersection point for one of the less refractory 
condensates (such as forsterite), and exponential decays 
on both sides (si = 6; S2 = 10). si = 6 is steep, but 
comports with the expected thinness of single clouds of 
refractory condensates having large particle radii. Mod- 
els E and D are similar. 

Our default value of toq for the composite cloud model 
E is where Tg for the most refractory condensates (the 
calcium aluminates) equals the atmospheric tempera- 
ture. For simplicity, in our fiducial model E we have 
set To equal to a constant 2300 K (Fig. Our toi for 
all models is taken to coincide with the intercept with 
the forsterite condensation curve (the forsterite TOc, at a 
temperature between 1500 and 1900 K; see Fig. 

Hence, the remaining parameters of our cloud mod- 
els are the modal particle radius (ag), and the optical 
properties of the representative condensate. We have 
studied the consequences of using enstatite, corundum, 
and forsterite for the latter and have settled on forsterite 
(Mg2Si04) for this generation of dwarf models. Its abun- 
dance should be one of the largest in the mix, and, as we 
will see ([Q, when the particle radii are near ~50— 100 
fim , the spectral fits over most of the relevant Tgff and 
gravity range are rather good. 

4. EXPLORATION OF PARAMETER DEPENDENCES 

To clarify the role of cloud extent as described in H3.ll 
Fig. [^compares theoretical spectra for cloud distribution 
models A, B, C, D, and E at TcffS of 1100 and 1700 K, a 
gravity (g) of 10^ cm s"^ and a modal particle radius of 
30 /xm . This figure does not represent our preferred set 
of models, but rather is meant to portray the generic dif- 
ferences between the cloud models described in i l3.1l At 
Tcff — 1100 K, models C, D, and E are very similar, since 
the cloud base where these models differ is all but invis- 
ible for this Teff . These models also maintain the large 
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band peaks. Models A and B, which allow the cloud to 
extend to the lowest pressures, are severely affected by 
the cloud opacity, and result in very flat spectra and a 
redistribution of flux from short to long wavelengths that 
is completely inconsistent with observations (see ^Sj) . On 
the bottom panel of Fig. 13 we show the corresponding 
models for Teff = 1700 K. Here, the near degeneracy be- 
tween models C and D/E is broken and for model C we 
can discern in the stronger J and H peaks the absence of 
cloud below the forsterite condensation curve. Clearly, 
the extent of the cloud(s) does make a difference. 

In Fig. 121 we compare models A, B, C, D, and E for 
30-/im particles on an M j versus J — K color-magnitude 
(HR) diagram, but for a wider range of Teg s than pro- 
vided in Fig. [5] Included are both the cloud-free (clear) 
curve and a collection of measured M, L, and T dwarfs 
with parallaxes. All the cloud models (except model A) 
reproduce the general behavior of the M— >L transition 
and the C, D, and E models comfortably inhabit the 
early-to-mid L dwarf regime. The presence of clouds 
clearly defines the Ls. In addition, models C, D, and 
E are degenerate below Teff ^^1400 K at this gravity and 
reach the measured mid-T dwarf region to the blue in rea- 
sonable fashion. Note that the extended cloud models A 
and B never make the transition to the "clear" behavior 
of the late T dwarfs and they are off in J — ivT by as much 
as 3 magnitudes at low Toff . However, though as Fig. 
12 demonstrates the corresponding spectra for models A 
and B are disfavored, even these models transverse the 
L-type spectral band in reasonable fashion. The same 
can be said for the models of Chabrier et al. (2000) (to 
which models A and B are closest in their particulate 
spatial distribution) for Mjs from ~11 to ~14.5. How- 
ever, depending upon the isochrone chosen, the Chabrier 
et al. models later than ^^15.0 magnitudes are redder 
than models A and B by as much as '--^0.5 to ~2.0 mag- 
nitudes. This is in part a consequence of their use of an 
ISM particle size distribution. 

Furthermore, and importantly, as shown in Fig. 13 the 
observations of early- to mid-T dwarfs indicate a substan- 
tial brightening in J, peaking near T3. If in fact this ob- 
servational trend is denoting brightening along an evolu- 
tionary track or along a near-constant gravity/metallicity 
trajectory, and not an incompletely sampled scatter in 
the known population of brown dwarfs with parallaxes 
or the effects of binarity, then models C, D, and E, 
though they handle the early-to-mid Ls and the mid- 
to-late Ts reasonably well, do not adequately reproduce 
the L^T transition region. The diminution of the effects 
of cloud optical depth with decreasing Toff is occurring 
more rapidly in the real objects than in the models. We 
will discuss this more in S^and [IHl 

Figure 0| portrays the dependence on the exponential 
cutoff parameter, Si, of the theoretical spectrum of an L 
dwarf with a Tcs of 1500 K, a gravity of 10^ cm s~^, a 
modal particle radius of 30 /im , and the model E par- 
ticle density distribution. The purpose of this diagram 
is to demonstrate that the upper (low-pressure) particle 
density scale height, a quantity that is very loosely con- 
strained by first principles and that may be chosen by 
modelers with little physical motivation, makes a differ- 
ence of as much as ^20% in the near infrared fluxes. Ig- 
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the si the thinner the upper cloud deck. For our sub- 
sequent model E spectra, we chose si = 6 to reflect our 
prejudice that the refractory clouds made up of large par- 
ticles settle more compactly than terrestrial water clouds 
(Ackerman & Marley 2001). 

Next, we portray in Fig. [3 an HR diagram similar 
to Fig. 121 but as a function of modal grain particle 
size/radius from 3 to 100 microns. For comparision, the 
cloud-free trajectory is included. Toff s from 700 K to 
'--^2000 K are represented (see caption) and model E with 
si = 6 is employed. A constant gravity of 10^ cm s~^ and 
forsterite optical properties are assumed. Data for M, L, 
and T dwarfs with known parallaxes are superposed. 

Figure El shows that, at least in Mj/(J — K) space, 
particle size makes little difference for the early- to mid- 
L dwarfs, but a very large difference in the turnaround 
elbow and in the Mjs for the early- to mid-T dwarfs. 
The trajectory with a larger particle size approaches the 
anomalous T3 brightness, but fails to include the lat- 
est, reddest L dwarfs. A small particle size can "fit" 
the latest L dwarfs, but completely misses the early T 
dwarfs. In short, all else being equal no one particle size 
and gravity combination can be made to fit all the data, 
even from the perspective of this lone HR diagram. One 
might conclude from Fig. [3 that a systematic increase 
in oq with decreasing Teff might be indicated (Knapp et 
al. 2004), but how this is accomplished by nature would 
remain to be determined. In fact, the simple convective 
cloud models of Cooper et al. (2003) and Ackerman & 
Marley (2001) predict that particle size slowly decreases 
with decreasing Toff . However, a fast transition to large 
particles as Teff decreases, or at some critical Toff , due 
perhaps to the onset of runaway particle growth and pre- 
cipitation in convective flows with large shears (Shaw 
2003; Kostinski & Shaw 2005), would be worth investi- 
gating. 

However, as is demonstrated by Fig. El for a set of 
Toff = 1500 K and g = 10^ cm s~^ models as a function of 
modal particle radii and by a basic knowledge of L dwarf 
spectra (as well as of the figures in , small particle sizes 
below ~30 microns result in theoretical spectra that are 
much too smooth, that don't provide the high contrasts 
across the water bands generically observed for objects 
later than mid-L dwarfs, that suppress the K I doublet 
in the J band, and that ruin the CO and methane fits 
between 2.2 and 2.4 microns. Hence, small values for ao 
are strongly disfavored for L and T dwarf models. The 
upshot is that while larger particles (> 40 yitm ) seem 
favored from general spectral considerations, they make 
fitting the reddest late L dwarfs problematic. From now 
on, we will take ag = 100 /imfor our model E (with 
si — 6) calculations as the best representative parameter 
combination and see where that can lead us. 

In picking a single particle size for our reference model, 
we are not only acknowledging, as we hope others would 
also do, our current ignorance of the physics and sys- 
tematics of clouds in brown dwarf atmospheres, but our 
observation that such large particle sizes, whatever the 
grain kinetics, are favored by the data. We find in fact, 
as the following sections show, that our baseline model 
fits the data rather well, and certainly no worse than any 
other model now in the literature. The reader is invited 
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limitations of our baseline model. 

A cloud model must provide the values and systematics 
with composition, Tcff , and gravity of the modal particle 
size, the particle size distribution, the spatial distribution 
(with altitude) of cloud properties, the optical constants 
of all the relevant grains as a function of wavelength 
(with the effects of grain mantling), the effects of grain 
asphericity, the patchiness of the clouds, the effects of 
convective overshoot, and the variations with longitude 
and latitude. Though there have been some good starts 
at framing the discussion of some of these components 
(Helhng et al. 2001, 2004; Woitke & Helling 2003,2004; 
Ackerman & Marley 2001; Cooper et al. 2003), we be- 
lieve these issues are still to be adequately addressed. 
Therefore, in this paper we provide our parameter study 
and a phenomenological reference model, avoiding what 
does not work, while being guided by what seems to, all 
without pretending to have achieved an understanding 
of cloud and dust physics. 

Note that the small abundances of the cal- 
cium/aluminium refractories formed in the radiative zone 
near the M— >L transition, coupled with the higher per 
particle opacities of the small particles expected in such 
zones, result in a total optical depth that is coincidently 
comparable to that for the 100- micron forsterite particles 
assumed in our baseline model. Hence, the total optical 
depth of large particles of high-abundance species (such 
as forsterite) can crudely mimic the total optical depth 
of the smaller/low-abundance particles expected in ra- 
diative zones. 

Before we transition to a discussion of our fiducial 
dwarf model set, we present in Fig. [7|the dependence 
on gravity of Toff trajectories in Mj versus J — K space 
(for flQ = 30, 100 /im). Each color represents a different 
gravity (either lO''-^, 10^ ", or 10^-^ cm s~^), with the 
larger gravities to the left. For each color, the line to the 
left is the larger particle size (100 /xm). 

The range of gravities employed for Fig. drepresents a 
reasonable range of dwarf masses, from ~15 A'/j to ~70 
Afj (Burrows et al. 1997). That this range of gravities 
nicely covers the data points in both the L and T dwarf 
regimes strongly suggests that the observed widths of 
the L and T dwarf bands are explained by gravity vari- 
ations, that g is the "second parameter." Though one is 
still unsure about the latest Ls and the brightening in J 
near T3, this conclusion seems reasonable and has also 
been reached by Burrows et al. (2002b) and Knapp et 
al. (2004) using different arguments. However, it is still 
possible that binarity can be a factor in the width of the 
L and T spectral bands in the HR diagrams. 

5. TEMPERATURE/PRESSURE PROFILES 

Now that we have explored various cloud parametriza- 
tions, and settled upon one of them to comprise our de- 
fault set (model E, with oq = 100 /im , si = 6, and 
forsterite's complex index of refraction versus wavelength 
behavior), let us return to a discussion of the regions and 
structure of the generic L or T dwarf atmosphere. To do 
this, we provide in Figs. |H1 and El temperature/pressure 
profiles for different TeffS spanning the L and T dwarf 
Tcff range, at gravities of 10^'^ cm s~^ and lO'"^ " cm 
s~^, respectively. The red regions are convective, the 
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of the average photosphere. The dashed lines indicate 
the approximate positions of the inner and outer cloud 
boundaries, as described in i i3.1l 

As Figs. ElandElshow, at TcffS from 2200 K to 1900 K 
clouds are radiative and the photospheres are in the ra- 
diative zone. At lower ToffS , from ^--^1800 K to ~1600 K, 
the precise value of which depends upon gravity (and, in 
fact, particle size and metallicity), an isolated outer con- 
vective zone emerges that coincides approximately with 
the position of the condensate cloud. The increasing 
opacity of the particulate matter trips local convection, 
but at the same time flattens the T/P profile below this 
zone at slightly greater pressures. The upshot is an iso- 
lated radiative zone between two convective zones. Even- 
tually, near ^--^1600 K for the default model depicted here, 
the two convective zones merge. What happens meteoro- 
logically and dynamically when they merge might make 
for an interesting study; we will return to this later. Af- 
ter the two convective zones merge, the outer convective 
boundary is always near 1500—1600 K and ranges in pres- 
sure from ^0.5 to ~20 bars for the lower gravity and from 
~1 to ^80 bars at the higher gravity. As Figs. HI and |^ 
demonstrate, the pressures at the cloud base for the two 
model sets range from ~0.3 to ~20 bars for g=l{f -^ cm 
s~^ and from ~0.7 to '^100 bars for g=10^ cm s~^ At 
a given atmospheric temperature, increasing gravity or 
decreasing metallicity increases the local pressure. Cen- 
tered at a gravity of 10^ cm s~^ and solar metallicity, 
at a given T a one order of magntiude range in either 
gravity or metallicity translates into a variation of about 
an order of magnitude in pressure. 

The fact that the pressures at the cloud base inex- 
orably rise with decreasing Toff implies that the optical 
depth of the cloud increases dramatically with decreasing 
Toff , all else being equal. The Rosseland optical depth 
of the condensate clouds in the models in Figs. |Hland|ni 
range from of order unity to one hundred. The clouds 
get thicker with later spectroscopic type from the early 
L dwarfs to the late T dwarfs. 

However, as the positions of the asterisks indicate, the 
photosphere ranges from '--^0.2 to ~0.6 bars at the lower 
gravity and from ~0.5 to ^^2.0 bars at the higher grav- 
ity. Its position is tightly constrained to lie in a narrow 
pressure range around a bar. As a result, it recedes from 
the cloud base and cloud tops with decreasing Toff ; as 
Figs. IHl and [HI demonstrate, the top of the convective 
zone is progressively buried and separates from the "pho- 
tosphere." This is the origin of the waning influence as 
Toff decreases of clouds in dwarf atmospheres and spec- 
tra. Note that as Toff decreases, the photosphere is first 
in a radiative zone (early L), then the convective zone 
(mid L), and finally in the outer radiative zone. Figure 
^1 depicts the Rosseland depth above the clouds due to 
function of Teff for various gravities and shows 
how as Toff decreases the cloud tops are more and more 
obscured by the cloud-free, gas-dominated regions. 

On the basis of the discussion above, one might con- 
clude that the L dwarfs, the transition from L to T, and 
the brightening in the J band depicted by the data plot- 
ted on Figs. 13 |S1 and [3 all naturally follow and can be 
explained by the phenomenology represented in Figs. |S1 
andEl and most of this is true. However, the L— >T tran- 
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from the photosphere, or the particle size increases, faster 
than is imphed by the models of Figs. |H1 and |5| what- 
ever their other strengths. It has been suggested that 
the clouds break up and spawn holes that widen with 
decreasing Toft (Burgasser et al. 2004). Since the clouds 
inexorably thicken as Toff decreases, we do not favor this 
break up hypothesis. However, given our current igno- 
rance of cloud meteorology, dynamics, and grain growth, 
not the least at the cloud tops that for late-L and early- 
T dwarfs are closest to the photospheres, we can not 
at present eliminate this scenario. Be that as it may, 
we can compare the spectral and color models we de- 
rive with representive observed dwarf spectra to conclude 
that apart from the issue of the L-^T transition, our new 
models reproduce the observations rather well. Next, 
we demonstrate this, using only the coarse model set in 
Toff and gravity space we have generated for this paper. 

6. SPECTRAL COMPARISONS 

Figures [TTlll2l andll^lsummarize the spectra in the op- 
tical and near infrared associated with our baseline model 
set as a function of Toff (Fig. ^] at solar metallicity) , 
gravity (Fig. 1121 for our three gravities), and metallicity 
(Fig. ^1 for three metallicities)^. These models fully in- 
corporate an algorithm for condensation clouds and the 
latest gas-phase opacities. As Fig. indicates, and 
as one would generally expect, lower metallicity models 
are bluer (compare the green, red, and blue curves at 
the Y jZ^ J, and H peaks). However, some of the non- 
monotonic behavior in, for instance, the K band can be 
explained by the fact that the CO/CH4 ratio, the atmo- 
spheric pressures, and the contribution of clouds are all 
metallicity-dependent, sometimes in countervailing ways. 

In addition to incorporating the evolving effects of 
grains as Teg decreases, these spectral models track the 
transition from CO to CH4, particularly manifest in the 
K band in the 2.3 /imto 2.2 /im switch, the deepening 
of the methane feature in the iJ band, the appearance 
(then disappearance) of FeH and CrH features near ~1.0 
^mand ~0.85 /xm, the emergence of the neutral Na, K, 
Cs, and Rb alkali features below 0.9 /xm, the emergence 
of the Y I Z band peak (near 1.05 /xm), the disappearance 
of TiO and VO, and the resculpting of the J, if, and K 
bands. The shape change in the K band due in part to 
collision-induced H2 absorption is particularly striking. 
Also seen is the evolution of the K I features in the J 
band near 1.25 /imand near 1.17 \iva.. 

Without attempting detailed fits, we compare in Figs, 
rm and El the closest representatives of our model set 
at solar metallicity with four observed spectra for the 
LI dwarf 2MASS J03454316-h2540233 (Kirkpatrick et 
al. 1997; Leggett et al. 2001), the L5 dwarf SDSS 
J05395199-0059020 (Leggett et al. 2000; Geballe et al. 
2002), the T4.5 dwarf SDSS J02074248+0000562 (Tsve- 
tanov et al. 2000; Geballe et al. 2002), and the T8 dwarf 
Gliese 570D (Burgasser et al. 2000; Geballe et al. 2001). 
All the measured spectra are for objects with measured 
parallaxes (Vrba et al. 2004) and Figs. El and El are 
absolute, not relative, comparisons. Had we generated a 
new set of models with a finer gravity. Tog , and metal- 
licity grid, we could have obtained even better fits, but 



have left that exercise to a future paper. 

For the LI dwarf 2MASS J0345, though the H-hund 
comparison is slightly problematic, the correspondence 
between theory and observation in the shape of the spec- 
trum from 1.8 to 2.3 /imand in the J band is quite 
good. For the L5 dwarf SDSS J0539, the correspondence 
throughout the spectrum depicted is excellent, and even 
the FeH Wing- Ford feature at ~0.99 /imis reasonably 
well reproduced. Except for shape discrepancies in the 
H and K bands, the correspondence between theory and 
observation from 0.9 /imto 1.6 /xmfor the T4.5 SDSS 
J0207 is striking. Furthermore, for the T8 dwarf Gliese 
570D, except for deviations at the Y/Z peak and in the 
H band (the latter due, again, to problems with the ex- 
tant methane database), the "fit" with theory is remark- 
able. Note that all these comparisons are being made for 
spectra plotted linearly, not in the logarithm, and are, 
therefore, that much more encouraging. 

The successes for the L dwarf comparisons seen in Fig. 
Elsuggest that for the early to mid-L dwarfs our conden- 
sation cloud model is rather good, surprisingly so given 
the possible variations on cloud modeling. Furthermore, 
Fig. 1151 suggests that for the mid- to late T dwarfs the 
models are also good. For all four representative dwarfs 
depicted, the ToffS and gravities shown on the plots are 
useful first guesses for these quantities. As Marley, Gush- 
ing, & Saumon (2004) suggest, for T dwarfs later than 
about T5 we find that the clear models without any 
direct effects of clouds fit the observations in the near 
to mid-infrared slightly better than the cloudy models. 
Marley, Gushing, & Saumon (2004) claim that this de- 
marcation is near Teff s of 1200 K. However, we determine 
that the effects of clouds, particularly from the optical to 
~1.1 microns (Burrows et al. 2002), completely disap- 
pear at ToffS cooler than ^1100 K. This slight difference 
highlights the problems that attend the L^T transition. 
Beyond TeffS of 1100 K we discern little difference at any 
wavelength for which comparisions can be made between 
the Marley, Gushing, & Saumon (2004) models and the 
models presented in this paper. 

7. COLOR-MAGNITUDE RESULTS 

If there were to be a brightening due to the diminu- 
tion of the effects of clouds, it is not surprising that this 
would be most manifest in the J and Y/Z bands. Figure 
1161 provides a comparison between the forsterite and the 
solar-metallicity gas opacities for two particle sizes and 
a typical atmospheric thermodynamic point. The gas- 
phase opacities show the classical variation due to the 
water absorption troughs that define the near infrared 
photometric bands. At ^1.05 /xmand ~1.25 /im, the 
grain opacities are most competitive with the gas, indi- 
cating clearly that it is in the Y/Z and J bands that 
the recession to depth, thinning, or disappearance of the 
clouds would first be obvious. Hence, it is all the more 
perplexing that our default models can not reproduce 
the ~1.0 magnitude rise seen in Mj versus J — K near 
subtype T3, if "rise" it is. Nevertheless, it is instructive 
to see just what the new models do predict. Figure ITTl 
is a Mj versus J — K HR diagram with Tog trajectories 
for all our fiducial models for the three gravities (lO'*'^, 
10^ °, and 10^-^ cm s"^) and three metallicities ([Fe/H] 
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also included on the plot. A given color denotes a given 
metallicity, the leftmost models being the most metal- 
poor. The green curves are for solar metallicity. For a 
given triplet, the leftmost model is the highest gravity. 

Consonant with the discussion in 21 the early to mid- 
L dwarfs and the mid- to late-T dwarfs fit quite well. 
The thicknesses of both the T and L dwarf bands are 
well recreated by variation in gravity, with perhaps a 
slight metallicity dependence. The gravity dependence 
in the L— s-T transition is expected to be weak, while the 
metallicity dependence is strong, as large as could be 
produced by a variation in particle size of a factor of ten. 
The latest L dwarfs are still outliers in the theory, for 
all TcffS , gravities, and metallicities. It is indeed ironic 
that to put models among these latest L dwarfs requires 
higher cloud opacities, while to explain the anomaly near 
T3 requires lower cloud opacities, emphasizing that the 
swing from L9 to T3 is very rapid, even more rapid that 
the discussion and plots in S|5l provide for. 

As Fig. El indicates, we could fit the J "brightening" 
around T3 with lower metallicities, only a factor of ^3 
below solar. However, as Fig. 1171 also indicates, such a 
metallicity, though it might help with some T dwarf fits, 
is not consistent with most of them. Furthermore, if the 
observed trajectory from L to T is not very incomplete, 
then the metallicity explanation would be hard-pressed 
to explain the empty regions in Mj/(J — K) space that 
should be occupied by low- metallicity dwarfs. 

The corresponding diagram for versus J — K (Fig. 
I18|l is more encouraging, with the green solar-metallicity 
lines inhabiting the L dwarf, the L— >T transition, and the 
T dwarf regions rather well. However, the latest L dwarfs 
are still left out. Note that Figs. ITTI and [THl suggest that 
the low-metallicity ([Fe/H]<-1.0) subdwarf L dwarfs will 
swing to the blue earlier than solar-metallicity L dwarfs 
and, hence, that the extent of the subdwarf L dwarf fam- 
ily in color and Tes space will be smaller. At low enough 
metallicity, the L dwarfs might altogether disappear as a 
class. 

Figure El depicts the J — H versus H — K color-color 
diagram for a subset of the new models and explains 
a curious behavior. As before, each line is for a given 
gravity and all models are for solar metallicity. Data for 
T, L, and M dwarfs and a cloud-free trajectory (orange) 
are also shown. As the data indicate, there is a distinct 
clump of M dwarfs. To the red in both colors can be 
found the L dwarfs, also as a clump. However, the T 
dwarfs are found on the opposite side from the L dwarfs, 
re-traversing the M dwarf realm and extending to the 
blue. The clear models do not extend redward of the 
M dwarfs. However, the cloudy models do show this 
behavior, curling up to the L dwarfs, then back through 
or near the M dwarfs, and on to the blue T dwarfs. This 
sinus behavior is a consequence of the appearance, then 
disappearance of the clouds from the photosphere. The 
best fit is for models (the dashed curves) for which the 
CH4 opacity in the H band has been artifically increased 
by a factor of 2 to account for the missing CH4 hot bands. 
Our default model set does not incorporate this artifical 
enhancement, but we present such models here to show 
the better fits that are possible with the same opacity 
increase as is required for better spectral fits to measured 
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sus Toff trajectories for a variety of metallicities and grav- 
ities. The basic behaviors of the entire M, L, and T dwarf 
families are encapsulated here, though the slope from the 
late L dwarf bump to the T dwarfs at lower Tcffsis too 
shallow to explain the anomalous brightening in J. 

Finally, in Fig. |2] we portray model trajectories for 
the various metallicities and gravities of our new model 
set in the M^/ versus i'-z' color-magnitude diagram^. L 
and T dwarf data from Knapp et al. (2004) and for which 
there are parallaxes are superposed. The theory makes 
a passable job at following the dwarf points, though for 
the few T dwarfs included a better fit is indeed with 
the lower metallicity lines, while for the latter L dwarfs 
a slightly higher metallicity would be preferred. More 
data in these bands for dwarfs with parallaxes is clearly 
needed. 

Figures 1111 through 1211 represent a new generation of 
models and the degree to which they reproduce current 
L and T dwarf data. All in all, the new models are an 
advance and they can be used with profit to extract phys- 
ical information from observed spectra, particularly for 
the early- to mid-L dwarfs and the mid- to late T dwarfs. 
However, a number of questions remain unresolved, in 
particular the character of the J band brightening and 
reason for the redness and dimness of the latest L dwarfs. 

8. CONCLUSIONS AND DISCUSSION 

Using an updated gaseous opacity database and a 
model for clouds in L and T dwarf atmospheres, we have 
generated and presented a new series of spectral models 
spanning the Tcfi range from 2200 K to 700 K, as a func- 
tion of gravity and metallicity. The correspondence with 
observed spectra and infrared colors for early- and mid-L 
dwarfs and for mid- to late-T dwarfs is good. We find 
that the width of both the T and L dwarf branches is nat- 
urally explained by reasonable variations in gravity and, 
therefore, that gravity is indeed the "second parameter" 
of the L/T dwarf sequence. Toff being the first. We have 
explored the dependence of theoretical dwarf spectra and 
HR diagrams upon various cloud parametrizations, such 
as particle size, cloud scale heights, and spatial distribu- 
tion. The results are a cautionary tale that our current 
ignorance of detailed cloud meteorology renders ambigu- 
ous the extraction of various physical quantities such as 
Toff and gravity for mid-L to early-T dwarfs. We es- 
timate that errors of '^50—100 K in Toff and '^0.3 in 
logj^Q (7 are likely, particularly near Toff ~1100— 1500 K. 
Nevertheless, for decreasing Toff , we capture with some 
accuracy the transition from CO to CH4 in the K band, 
the FeH and CrH features at ~1.0 //mand ^0.85 /im, 
the emergence of the neutral Na, K, Cs, and Rb alkali 
features, the growth of the Y /Z band peak, the disap- 
pearance of TiO and VO in L dwarfs, the transforma- 
tions of the J, H, and K bands, and the evolution of 
the K I features near 1.25 /j,mand near 1.17 /im. More- 
over, we speculate that the subdwarf branch of the L 
dwarfs will be narrow, and, that for low enough metal- 
licity, the L dwarfs will disappear altogether as a spectro- 
scopic class. Furthermore, we note that the new, lower 
solar oxygen abundances of AUende-Prieto, Lambert, & 
Asplund (2002) produce better fits to brown dwarf data 
than do the older Anders & Grevesse (1989) values. How- 
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ever, we do not reproduce the "brightening" in J, nor the 
dimness of the latest L dwarfs. 

What could be the reason for the problems near T3 
and L8/9? The J and YjZ band brightening will re- 
quire a more rapid decrease in atmospheric cloud opac- 
ity than occurs naturally in all the uniform cloud models 
published to date, including our own. Those of Tsuji, 
Nakajima, & Yanagisawa (2004) (their "UCM" set) do 
not reproduce either anomaly - neither do the models 
of Ackerman & Marley (2001), Marley et al. (2002), 
Chabrier et al. (2000), or AUard et al. (2001). Burgasser 
et al. (2004) have postulated the break up of the clouds 
near a T^ff of 1200—1300 K and the appearance of holes, 
whose filling fraction increases across the L— >T transi- 
tion until that fraction is unity. A virtue of this model 
is the natural explanation of the apparent resurgence of 
the FeH features in the early- to mid-T dwarfs (Bur- 
gasser et al. 2004; McLean et al. 2003; Gushing, Rayner, 
& Vacca 2005). The FeH abundances near the photo- 
spheres should be waning; holes could allow us to see 
more deeply to the higher-temperature regions in which 
the FeH abundance is large. The break-up model might 
imply larger temporal variations in the FeH, J, and YjZ 
band fluxes for early-T and late-L dwarfs, and such varia- 
tions can be looked for. However, no one has yet provided 
a plausible physical explanation for this break-up behav- 
ior, and the inexorable thickening of the cloud material 
with decreasing Toff would seem (superficially?) to miti- 
gate against this. Curiously, the interplay between rapid 
rotation and silicate clouds for some critical parameter 
regimes might result in banding, such as is seen in the 
atmosphere of Jupiter, and this may be playing a role 
in the observed J-band brightening. However, how such 
a cloud-rotation coupling could affect a transition from 
profound cloud opacity to full clarity is not immediately 
obvious. Be that as it may, the problem of fitting the 
latest L dwarfs on the infrared HR diagrams, while si- 
multaneously also fitting their spectra, is not explained 
at all by the filling factor hypothesis''. In fact, it might 
exacerbate it. 

It could be that the anomalous brightenings in J 
and YjZ are the result of crypto-binarity. The binary 
fraction of T dwarfs is not negligible and ~0.75 mag. 
(2.5 X logj^Q 2) is near the magnitude of the few excesses 
measured. It could also be, as we show in Figs. El that 
low metallicity is the culprit. However, both explana- 
tions seem rather ad hoc and do not explain the gaps 
and empty regions in the Mj/(J — HR diagram cur- 
rently inferred. The redness and low J-band fluxes seen 
for the latest L dwarfs might be explained by errors in 
the parallaxes or photometry. These objects might also 
be explained by supersolar metallicities. However, both 
these explanations seem unlikely, particularly given the 
largish number of L dwarfs that now inhabit that region 
of phase space. 

We have shown that when calcium-aluminate, silicate, 
and Fe clouds first form they do so in the radiative re- 

Note that we can fit the J vs. J — K points of the late Ls with 
small modal particle sizes, but can not simultaneously reproduce 
their near-infrared spectra from 1.0 to 3.0 microns. Small particle 



gion, that as Toff decreases an isolated convective zone 
emerges, and that for even lower TeffS the two convec- 
tive zones join. It remains to be seen what happens to 
the particle sizes and cloud morphology when these re- 
gions merge. Shaw (2003) and Kostinski & Shaw (2005) 
investigate the dependence of runaway droplet growth 
and rainout on the presence in convective clouds of large 
velocity shears and on intermittency in the turbulence. 
Could the merger of the outer convective cloud with 
the inner convective zone lead to regions of such large 
shears, in which particle growth on the timescales avail- 
able is more rapid, and, hence, to very large particles? 
Could the merger lead to the irreversible partial flushing 
of cloud material into the interior? After the joining of 
the convective zones, is the timescale for grain growth too 
long for the convecting feedstock to avoid being dragged 
into the hot interior before forming opaque grains? Is a 
critical TcS /gravity threshold for rapid grain coalescence 
and growth reached, beyond which the average particle 
is too large to contribute significant opacity (Liu, Daum, 
& McGraw 2005)? Or does the scale height of the sili- 
cate cloud collapse at some Toff threshold? The answers 
to these questions require a multi-dimensional approach 
both to grain kinetics and growth and to convective cloud 
structures and motions, all properly coupled. 

Since the discovery of Gleise 229B ten years ago, the- 
orists have made great strides in modeling brown dwarf 
atmospheres and spectra. The thermochemical and spec- 
troscopic data for the relevant molecules have improved 
considerably, as have the codes needed to handle strat- 
ified atmospheres. Observers have identified whole new 
spectroscopic classes, and the observations have got- 
ten spectacularly better across the electromagnetic spec- 
trum. Despite this palpable progress, much remains 
to be done. Foremost are the creation of physically- 
motivated cloud models, the expansion of molecular line 
lists that even now remain incomplete, the discovery of 
what lies beyond the T dwarfs at still lower ToffS , the 
further calibration of evolutionary and spectral models 
with well-calibrated spectra, and the discovery and char- 
acterization of more eclipsing systems. Robust gravity 
and metallicity diagnostics need to be determined and 
pursued and the binary statistics of substellar-mass ob- 
jects need to be put on a firm foundation. We look for- 
ward to participating in the next generation of progress 
on these problematic, yet fascinating, objects that have 
only recently emerged from nowhere to reawaken stellar 
astronomy. 

The authors thank Mike Gushing, Adam Burgasser, 
Bill Hubbard, Jonathan Lunine, Christopher Sharp, Jim 
Liebert, Drew Milsom, Curtis Cooper, and Jonathan 
Fortney for fruitful conversations and/or help during 
the course of this work, as well as NASA for its finan- 
cial support via grants NAG5-10760, NNG04GL22G, and 
NNG05GG05G. 

sizes (< 20 livci) severely flatten late L dwarf spectra in ways not 
consistent with the observations. 
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Fig. 1. — Solar-metallicity condensation curves (temperature in Kelvin versus pressure in atmosplieres) for many of the most important 
refractory species thought to appear in the atmospheres of brown dwarfs and L dwarfs. The most refractory compound is the calcium 
aluminate grossite (CaAUOy = CaO + 2(Al203)). Corundum (AI2O3), as such, does not generally form. CaMgSi206 is diopside, Mg2Si04 
is forsterite, MgSiOa is enstatite, MgAl204 is spinel, and Ca2MgSi207 is akermanite. The dotted curves correspond to the refractory 
titanium compounds. Contrary to popular belief, condensed titanium is rarely in the form of perovskite (CaTiOs). The golden dashed 
curve is for condensed VO. Liquid Fe is the solid peach curve at a slightly shallower slope than those for the calcium/aluminum/magnesium 
condensates (solid colors). Included for comparison are the condensation curves for water (H2O, red dashed) and ammonia (NH3, green 
dashed). Notice how the condensation curves of the refractory compounds densely inhabit a narrow range of T/P space and that there is 
a noticably wide gap between this refractory band and water. See text for a discussion of the salient features of this figure. 
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Fig. 2. — top panel: Near-IR spectra of brown dwarfs with Teff =1100 K and g=10^ cm s~'^ for differing cloud shape prescriptions 
(A through E). A forsterite modal cloud particle size of 30 firais assumed. The spectra produced using cloud shape models C, D, and E 
(red, dashed green, and blue curves, respectively) result in very similar spectra, which largely overlap on this plot. In contrast, cloud shape 
models A and B produce much redder spectra (black and purple curves). For comparison, a cloud- free spectrum for the same effective 
temperature is denoted by the dashed orange curve, bottom panel: Similar to top panel, but for Toft =1700 K. At this higher Tcft , more 
spectral variation due to cloud shape is exhibited. In particular, model C is closer to the cloud-free case, due to the smaller cloud opacity 
relative to models D and E. 



14 




Fig. 3. — Absolute J magnitude versus J — K color for brown dwarf models with g=10^ cm s~^ for 5 different cloud shape prescriptions 
and one cloud-free case. The points along each model curve are in Toff intervals of 100 K, ranging from 700 K to 1400 K for Case A, 700 
K to 2100 K for Cases B, C, and D, and 700 K to 2200 K for Case E and the cloud-free case. Shown for comparison are L and T dwarf 
data from Knapp et al. (2004) and M dwarf data from Leggctt (1992). A forsterite modal cloud particle size of 30 microns is assumed. 
Even in an approximate sense, the cloud-free model set and model sets A and B do not come close to fitting the complete set of brown 
dwarf data in this color-magnitude space. On the other hand, model sets C, D, and E come closer to reproducing the trends seen in the 
data, although the brightening in J around the early T dwarfs is not predicted. The strong turn toward the blue in J — K with decreasing 
Tj,ff near 1400-1500 K is due to the natural deepening of the cloud position in these model atmospheres. Throughout the paper, we use 
the MKO photometric system for the J, H, and K bands (Tokunaga & Simons 2002; Tokunaga, Simons, & Vacca 2003). 
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Fig. 4. — Dependence of a model dwarf spectrum on the vertical extent of the condensate cloud. A sharper cloudtop cutoff (larger 
exponential cutoff parameter, si) results in a bluer infrared spectrum. Our cloud particle distribution model E is assumed. 
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Fig. 5. — Absolute J magnitude versus J — K color for brown dwarf models with g=10® cm s ^ . Four model sets, with modal forsterite 
particle sizes of 3, 10, 30, and 100 fira , along with an additional cloud-free set, are plotted. The points along each model curve are in 
Tcff intervals of 100 K, ranging from 700 K to 1900 K for the 3- and 10-/im models, 700 K to 2200 K for the 30- and lOO-^tm models, and 
700 K to 1500 K for the cloud-free models. Shown for comparison are L and T dwarf data from Knapp et al. (2004) and M dwarf data 
from Leggett (1992). At a given Tcff , model atmospheres with larger particles result in stronger J-band fluxes. The strong turn toward 
the blue in J — ii" with decreasing Tcff near 1400-1500 K is due to the natural deepening of the cloud position in these model atmospheres. 
The model E cloud shape function is employed. 
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Fig. 6. — Dependence of a model dwarf spectrum on the modal cloud particle size of forsterite. Relative to a cloud-free model (dashed 
orange curve), near-IR spectral peaks and troughs are reduced due to the presence of clouds. Smaller particles result in redder, smoother 
spectra. Cloud model E parameters have been assumed. 
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Fig. 7. — Dependence on surface gravity and modal particle size for models in Mj versus J — K color-magnitude space. Cloudy (model 
E) and cloud-free model sets are plotted. For each gravity, the larger modal particle size (100 fim) falls to the left of the smaller modal 
size (30 fim) in this color-magnitude space. At a given Tj,ff , a higher gravity leads to a bluer J — K color and a reduced J-band flux. The 
points along each model curve are in T^ft intervals of 100 K, ranging from 700 K to 1500 K for the cloud-free models, 900 K to 2100 K for 
the 30-AJm , g=10''-^ cm s-^ models, 900 K to 2000 K for the lOO-^m , g=W' '^ cm s'^ models, 700 K to 2200 K for both cloudy g=10'^-° 
cm s-2 models, 900 K to 2000 K for the SO-^m , g=10^-^ cm s-^ models, and 800 K to 2000 K for the lOO-^tm, g=10^-^ cm s'^ models. 
Shown for comparison are L and T dwarf data from Knapp et al. (2004) and M dwarf data from Leggett (1992). 
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Fig. 8. — Temperature-pressure profiles for a set of models with constant surface gravity, logjQ <; = 4.5 and particle size, 100 microns, 
for different effective temperatures ranging from T^g = 2000 K (leftmost curve) to 900 K (rightmost curve) . The position where the local 
temperature is equal to the effective temperature (which indicates an approximate location of the region of formation of the emergent 
radiation) is shown as an asterisk. The cloud bases are depicted as dashed lines; the convection zones are drawn in red; and the black 
dots show the position of the boundaries of the convection zone(s). Notice the occurence of two distinct convection zones for effective 
temperatures between 1600 K and 1800 K. 
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Fig. 9. — Temperature-pressure profiles for a set of models with constant surface gravity, log^Q (;(cms~^) = 5.0 and particle size, 100 
microns, for different effective temperatures ranging from T(,fi = 2200 K (leftmost curve) to 700 K (rightmost curve). The position where 
the local temperature is equal to the effective temperature (which indicates an approximate location of the region of formation of the 
emergent radiation) is shown as an asterisk. The cloud bases are depicted as dashed lines; the convection zones are drawn in red; and the 
black dots show the position of the boundaries of the convection zone(s). Notice the occurence of two distinct convection zones for T^ff = 
1700 K and 1800 K. 
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Fig. 10. — Rosseland mean optical depth of gas above cloud versus Tgff as a function of surface gravity. With decreasing effective 
temperature, the clouds become buried beneath an increasing column of gaseous opacity and so the infrared spectra and colors of these 
brown dwarfs approach those of cloud-free atmospheres. 
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Fig. 11. — A selection of model E spectra (assuming a modal forsterite particle size of 100 microns) normalized to 10 pc. Spectra for 
eight different effective temperatures from 700 K to 2100 K are shown, encompassing the range from late T dwarfs to early L dwarfs. A 
surface gravity of 10^ cm s~^ is assumed and the theoretical radii of Burrows et al. (1997) are used. 
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Fig. 12. — Model spectra at effective temperatures of 1700 K (thick curves), 1400 K (medium curves), and 1000 K (thin curves) for 3 
different surface gravities, 10^'^, 10^, and 10^'^ cm s~^ (purple, cyan, and black curves, respectively.) The theoretical brown dwarf radii of 
Burrows ct al. (1997) are used, and our model spectra are normalized to 10 pc. 
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Fig. 13. — Model spectra at effective temperatures of 1700 K (thick curves), 1400 K (medium curves), and 1000 K (thin curves) for 3 
different metaUicities, [Fe/H]=-0.5, 0, and 0.5 (blue, green, and red curves, respectively.) A gravity of 10^ cm s~^ was employed. The 
theoretical brown dwarf radii of Burrows et al. (1997) are used, and our model spectra are normalized to 10 pc. 
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Fig. 14.— Spectrum of 2MASS J03454316+2540233 (Leggett et al. 2001; Kirkpatrick et al. 1997), an LI dwarf, and SDSS J05395199- 
0059020 (Geballe et al. 2002; Leggett et al. 2000), an L5 dwarf. Both spectra have been normalized to a their absolute fluxes at 10 pc. 
Shown for comparison are model spectra, also normalized to 10 pc, from our coarse model grid, and using the theoretical radii of Burrows 
et al. (1997). These are by no means model spectral fits to the data. Rather, they are absolute comparisons between observational and 
theoretical results. Model E parameters are assumed. 
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Fig. 15.— Spectrum of SDSS J02074248+0000562 (Geballe et al. 2002; Tsvetanov et al. 2000), a T4.5 dwarf, and Gl 570D (Geballe 
et al. 2001), a T8 dwarf. As in Fig. 1141 the spectra liave been normalized to 10 pc for an absolute comparison with theoretical model 
spectra. 
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Fig. 16. — Abundance-weighted comparison of forsterite opa<;ity for 30-;imand 100-/im modal particle sizes with that of the total gas 
opacity at a pressure of 1 bar and temperature of 1500 K. In the Y jZ and J bands, forsterite can be a dominant opacity source, depending 
upon the depth of the cloud layer in the atmosphere. 
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Fig. 17. — Metallicity and gravity dependence of brown dwarf models in absolute J magnitude versus J — K color space for model E 
parametrizations. Cloudy models (assuming a forsterite modal particle size of 100 fira) and cloud-free models (thin curves) are shown for 
surface gravities of 10*'^, 10^", and 10^'^ cm s~^ and for metallicities, [Fe/H] = -0.5, 0, and 0.5 (blue, green, and red curves, respectively). 
Each color-coded set of 6 models (3 cloudy and 3 cloud- free) contains the 3 different surface gravities, with the leftmost cloudy and cloud- 
free curves being the highest gravity. The points along each model curve are in Tog intervals of 100 K, ranging from 700 K to 1500 K for 
the solar metallicity cloud-free models, 800 K to 1500 K for the nonsolar metallicity cloud-free models, 900 K to 2000 K for the g=10*'^ cm 
s~'^ solar-metallicity cloudy models, 800 K to 1900 K for the 3=10^" cm s~^ supersolar-metallicity cloudy models, and 800 K to 2000 K 
for all other cloudy models. A difference of a factor of 10 in metallicity carries roughly a full magnitude difference in J — K color. Shown 
for comparison are L and T dwarf data from Knapp et al. (2004) and M dwarf data from Leggett (1992). 



29 



S 




J-K 



Fig. 18. — Same as Fig. 1171 except in absolute K magnitude versus J — K color space. 
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Fig. 19. — J — H versus H — K color-color diagram. Assuming a modal forsterite particle size of 100 ^m , models sets for surface gravities 
of 10**'5, lO^", and lO^ '^ cm s are plotted (blue, red, and green solid curves, respectively). In addition, we have plotted a cloud-free 
model set for g = 10®'" cm s~^ (orange curve). Because our H-band opacity database is missing the "hot bands" of methane, we estimate 
the effects of the expected additional opacity by artificially enhancing the methane opacity in the //-band by a factor of 2 (blue, red, and 
green dashed curves). The points along each model curve are in T^ff intervals of 100 K, ranging from 900 K to 2000 K for the fiducial 
g=104-5 cm s-2 model, 700 K to 2200 K for the fiducial g=105 cm s-^ model, 800 K to 2000 K for the fiducial g=10^-^ cm s-^ model, 
800 K to 2000 K for the enhanced g=10'''5 cm s-^ and g=10^-^ cm s'^ models, 700 K to 2100 K for the enhanced 9=10^ ° cm s-^ model, 
and 700 K to 2100 K for the cloud-free model. Shown for comparison are L and T dwarf data from Knapp et al. (2004) and M dwarf data 
from Leggett (1992). A model E cloud shape function is employed. 
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Fig. 20. — Metallicity and surface gravity dependence of model E dwarf models in J — K versus T^g space. The J — K color tends to 
peak at its reddest value in the 1400-1600 K T^fi range. J — K varies significantly with both metallicity and gravity, but it clearly turns 
bluer at cooler effective temperatures. The bluest colors are exhibited for low T(,ff , low metallicity, and high surface gravity. 
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Fig. 21. — Metallicity and surface gravity dependence of dwarf models in absolute z' magnitude versus i' — z' color space. Sloan filter 
functions and zero points were used. Cloudy models (model E) with a forsterite modal particle size of 100 /^m are shown for surface gravities 
of lO*'^, lO^", and 10^'^ cm s~^ (dotted, solid, and dashed curves, respectively), and for metallicities, [Fe/H] = -0.5, 0, and 0.5 (blue, 
green, and red curves, respectively). Shown for comparison are L and T dwarf data from Knapp et al. (2004). 



